Exciting half-integer charges in a quantum point contact 
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We study a voltage-driven quantum point contact (QPC) strongly coupled to a qubit. We predict 
pronounced observable features in the QPC current that can be interpreted in terms of half-integer 
charge transfers. Our analysis is based on the Keldysh generating functional approach and contains 
general results, valid for all coherent conductors. 



The quantum point contact^ has become a basic con- 
cept in the field of Quantum Transport owing to its sim- 
plicity. Its common experimental realization is a narrow 
constriction that connects two metallic reservoirs. An 
adequate theoretical description for this setup is a non- 
interacting one-dimensional electron gas interrupted by 
a potential barrier. The barrier is completely character- 
ized by its scattering matrix. This enables the scattering 
approach to Quantum Transport 2 . This allows one to 
describe the average current through the QPC, as well 
as fluctuations away from this average, in terms of single 
electrons passing through the constriction 3 -. The strength 
of the scattering approach is its ability to describe not 
only traditional realizations of a QPC, but all coherent 
conductors, including diffusive wires and tunneling junc- 
tions. 

Despite the correctness of the non-interacting electron 
description, truly many-body quantum correlations do 
exist and are observable in a QPC. They manifest them- 
selves in the full counting statistics of electron transfers 3 
and allow for detection of two-particle entanglement 4 
through the measurement of non-local current correla- 
tions. This suggests that the observation of many-body 
effects in a QPC crucially relies on a proper detection 
scheme. In this Letter, we give an example of how an ap- 
propriate detector uncovers such non-trivial many-body 
effects as half-integer charges. 

We probe the QPC with a charge qubit. Such a de- 
vice has already been realized using single and double 
quantum dots. Previously, the QPC has been used as a 
detector of the qubit stated. We propose a scheme in 
which the roles are reversed. Provided the qubit and 
QPC are coupled strongly, the switching between the 
qubit states is accompanied by severe Fermi-Sea shake- 
up in the QPC. The d.c. current in the QPC is sensitive 
to the ratio of the qubit switching rates and thereby pro- 
vides information about these severe shake-ups. 

Before analising the system in detail, the following 
qualitative conclusions can be drawn. The qubit owes 
its detection capabilities to the following fact: In order 
to be excited it has to absorb a quantum e of energy from 
the QPC. Here e is the qubit level splitting, a parameter 
that can be tuned easily in an experiment by means of 
a gate voltage. The QPC supplies the energy by trans- 
fering charge from the high voltage reservoir to the low 
voltage reservoir. The transfer of charge q allows qubit 



transitions for level splittings e < qV, V being the bias 
voltage applied. 

We can assume that successive switchings of the qubit 
between its states |1) and |2) are rare and uncorellated. 
The qubit dynamics are then characterized by the rates 
T2i to switch from state |1) to state 2) and Ti2 from 
\2) to |1). The stationary probability to find the qubit 
in state |2) is determined by detailed balance to be P2 — 
r2i/(ri2 + T2i). This probability can be observed exper- 
imentally by measuring the current in the QPC. The cur- 
rent displays random telegraph noise, switching between 
two values I\ and 12- These correspond to the qubit be- 
ing in the state |1) or \2) respectively. The d.c. current / 
gives the average over many switches and is thus related 
to the stationary probability by I — (1 — P2)h + P2h- 
The values of I\ , I2 and / are determined through mea- 
surement and P2 is inferred. 

When the QPC and qubit are weakly coupled^, a sin- 
gle electron is transferee^. This liberates at most energy 
eV, implying that the rate T2i is zero when e > eV and 
the rate Ti2 is zero when e < —eV. The resulting P2 
changes from 1 to upon increasing e within the interval 
— eV < e < eV. Cusps at e = ±eV signify that charge e 
is transferred. [See Fig. (2a)] 

Guided by our understanding of weak coupling we can 
speculate as follows about what happens at strong cou- 
pling. Apart from single electron transfers, we also ex- 
pect the coordinated transfers of groups of electrons. A 
group of n electrons can provide up to neV of energy 
to the qubit. Therefore, peculiarities in P2 should ap- 
pear at the corresponding level splittings e = ineV, 
n = 1, 2, 3, . . pi£ However, it is not apriori obvious that 
these peculiarities are pronounced enough to be observed. 
The reason is the decoherence of the qubit states induced 
by electrons passing through the QPC. The Fourier trans- 
form of the qubit transition rate acquires an exponential 
damping factor e~ w ^, W~ x being the decoherence time. 
This smoothes out peculiarities at the energy scale W . In 
the strong coupling regime, the decoherence time is esti- 
mated to be short, W ~ eV/ . As a result, it is not clear 
whether peculiarities at neV are the dominant feature at 
strong coupling. 

Therefore, strong coupling of the QPC and the qubit 
requires quantitative analysis. We have reduced the 
problem to the evaluation of a determinant of an infinite- 
dimensional Wiener-Hopf operator. We calculated the 



2 



determinant numerically and found that peculiarities at 
multiples of eV are minute. Their contribution to p 2 does 
not exceed 10 -4 and is seen only at logarithmic scale and 
at moderate couplings. Instead, far more prominent fea- 
tures occurs at e — \eV ■ General reasoning does not 
predict this. Straight-forward energy balance arguments 
force us to conclude that qubit switching is accompa- 
nied by the transfer of charge e/2 through the QPC. 
This frees up energy eV/2, stimulating qubit transitions 
when e < eV/2. In other words, the qubit switching ex- 
cites a half-integer charge and simultaneously detects it. 
Fractional charge is known to occur in strongly interact- 
ing many-electron system s 11 ! 12 ' 13 in equilibrium. In con- 
trast to this, the electrons in the QPC can be regarded 
non-interacting except during the short time the qubit is 
switching. Our system is also unusual in that the half- 
integer charge is only produced during qubit switching 
and is not present in the equilibrium state. 

Let us now turn to the details of our analysis. The 
system is illustrated in Fig. JT]). The Hamiltonian for the 
system is 



H = f+U, |1) (l\ + (U 2 +s) |2> (2|+ 7 (|1) (2| + |2) <1|) (1) 



order in the tunneling amplitude 7 it is given by 

,0 

r 2 i =2 7 2 Re / dr e leT 

J — OO 



x lim tr 

U)-^> — 00 



e iH2T e -iHi(T-t )p gii*i(T-to) 



(2) 



This is the usual Fermi Golden Rule. The Hamiltonians 
H\ and H 2 are given by Hk = T + Uk and represent 
QPC dynamics when the qubit is held fixed in state |fc). 
The trace is over QPC states, and po is the initial QPC 
density matrix. The evaluation of the integrand is a spe- 
cial case of a general problem in the extended Keldysh 
formalism^ 4 -. The task is to evaluate the trace of a density 
matrix after "bra's" have evolved with a time-dependent 
Hamiltonian H-(t) and "kets" with a different Hamilto- 
nian H + (t). 



tr 



<7-+ —i f°° dtH+(t) rr- i f°° dtH-(t) 



(3) 



We implemented the scattering approach to obtain the 
general formula 



A = tr In 



S_ (!-/)+«+/ 



— tr In s_ 



(4) 



The operator T represents the kinetic energy of QPC 
electrons. The operator Uk describes the potential bar- 
rier seen by QPC electrons when the qubit is in state 
h = 1, 2 and corresponds to a scattering matrix §k in 
the scattering approach. (We use a "check" to indicate 
a matrix in the space of transport channels.) QPC elec- 
trons do not interact directly with each other but rather 
with the qubit. This interaction is the only qubit relax- 
ation mechanism included in our model. We work in the 
limit 7 — > where the inelastic transition rates Pi2.2i 
between qubit states are small compared to the energies 
eV and e. In this case, the qubit switching events can be 
regarded as independent and incoherent. 

Now consider the qubit transition rate T2i. To lowest 




FIG. 1: A schematic picture of the system considered. It 
consists of a charge qubit coupled to a QPC. The shape of the 
QPC constriction, and hence its scattering matrix, depends 
on the state of the qubit. The QPC is biased at voltage V. 
A gate voltage controls the qubit level splitting e. There is a 
small tunneling rate 7 between qubit states. 



The operators s± and / have both continuous and dis- 
crete indices. The continuous indices refer to energy, or 
in the Fourier transformed representation, to time. The 
discrete indices refer to transport channel space. The 
operators s± = s±(t)6(t — t') are diagonal in time. The 
time-dependent scattering matrices s±(t) describe scat- 
tering by the Hamiltonians H±(t) at instant t. (It is the 
hall-mark of the scattering approach to express quanti- 
ties in terms of scattering matrices rather than Hamil- 
tonians.) The operator / = f(E)S(E — E') is diagonal 
in the energy representation. The matrix f(E) is diago- 
nal in channel space, representing the individual electron 
filling factors in the different channels. A full derivation 
of Eq. Q will be given elsewhere. It generalizes similar 
relations published i n 15 i 16 . 

In order to apply the general result to Eq. @, the 
time-dependent scattering matrices s±(t) are chosen as 

s+(t) =h + 6(t - T)6(-t)(s 2 - Si), (5) 
S_ =§!■ (6) 

The QPC scattering matrices Si(s2) with the qubit in 
the state 1 (2) are the most important parameters of our 
approach. 

Without a bias-voltage applied, the QPC-qubit setup 
exhibits the physics of the Anderson orthogonality 
catastrophe 1 ^. For the equilibrium QPC, the problem 
can be mapped 1 ^ onto the classic Fermi Edge singular- 
ity (FES) problemi&iQiS .. The authors o^ effectively 
computed A in equilibrium. Our setup is simpler than 
the generic FES problem since there is no tunneling from 
the qubit to the QPC. As a result, not all processes con- 
sidered in 1 ^ are relevant for our setup. We only need 
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the so-called closed loop contribution. The relevant part 
of the FES result for our setup is an anomalous power 

law T^2i(s) = 9(— (i^ - ) f° r the equilibrium rate. 
Here E CX1 . is an upper cutoff energy. The anomalous ex- 
ponent a is determined by the eigenvalues of s\si— as 



s 2 1 s 1 



1 

4?r 2 



Tr In (§fSi) 



The logarithm is defined on the 

branch (— tt,tt]. For a one or two channel point contact, 
< a < 1. 

We now give the details of our calculation for the rates 
out of equilibrium. From Eq. ([2]) and Eq. ([4]) it follows 
that r 2 i(e) oc |7 2 | f™^ dr e - l£T Det Q {v) {t). For positive 
times r, the operator Q^ v \t) is defined as^. 



Q (V) (r) = l 



(s 2 ^1 



i)n(r)/ 



(V) 



(7) 



while for negative t, Q^ v \t) = Q^(—t)^ The time- 
interval operator II(r) = 5(t — t')9(t)8(T — t) is diagonal 
in time and acts as the identity operator in channel space 
for times t = t' G [0, r] and as the zero-operator outside 
this time-interval. 

For the purpose of numerical calculation of the deter- 
minant we have to regularise Q^ v \t). This is done by 
multiplying with the inverse of the zero-bias operator to 

define a new operator Q(r) = Q {q) (t) _1 Q (v) {t). Its 
determinant is evaluated numerically. The rate r2i(e) 
at bias voltage V is then expressed as the convolution 
r 21 (e) = / %T e £(e - e')P{e') of the equilibrium rate 

and the Fourier transform of P{t) — Det Q^ v '(t), that 
contains all effects of the bias voltage V. 

We implemented this calculation numerically, and 
computed the probability P2 to find the qubit in state 
2). Details of our numerical method are presented in 
Appendix A. Our main results are presented in Fig. (2). 
We used 2x2 scattering matrices parametrized by 
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FIG. 2: The occupation probability p2 of qubit state j2). At 
weak coupling between the QPC and qubit, (Fig. a, b) the 
transfer of a single electron with charge e is detected. Pecu- 
liarities at ±eF/2 at strong coupling (Fig. c, d) constitute 
the detection of half integer charges e/2. Scattering matri- 
ces were parameterized as in Eq. [8] Fig. a, b, c and d 
respectively correspond to (f> = 7r/16, 7r/4, 7n/10 and 47r/5. 



cos q. 
i sin < 



i sin < 

cost? 



(8) 



and repeated the calculation for several <j) S [0, 7r]. Small 
(j) corresponds to weak coupling. The curve at 4> = 
is almost indistinguishable from the perturbative weak 
coupling limit discussed in the introduction. Cusps at 
±eV indicate that qubit switching is accompanied by the 
transfer of charge e in the QPC. 

The increasing decoherence smoothes the cusps for the 
curve at <fi = 7r/4 (2b). When the coupling is increased 
beyond <fi = ir/2 steps appear at ±eV/2 (c). This implies 
charge fractionalization e — > e/2. Further increase of the 
coupling results in a sharpening of the steps (d) . 

Known mechanisms of charge fractionalization do not 
seem to provide an immediate explanation of our find- 
ings. The Quantum Hall mechanismii does not give even 
fractions while the instanton mechanism^ requires a qua- 
siclassical boson field. There is an indirect analogy with 
the model of interacting particles on a ring threaded by 
a magnetic flujsi^. There, one expects that the energy 
eigenvalues are periodic in flux with period of one flux 
quantum. However, the exact Bethe-Ansatz solution 13 
reveals a double period of eigenvalues with adiabatically 
varying flux. This is a signature of half-integer charge 
quantization. 

For our non-equilibrium setup, energy eigenvalues are 
not particulary useful. The natural eigenvalues to de- 
scribe the phenomenon are those of the oprator Q( v '(t). 
They depend on the parameter eVr which is an analogue 
of flux. The product of the eigenvalues, i.e. the deter- 




1 2 3 4 5 6 1 

t [2Tr/eV] 

FIG. 3: The behavior of eigenvalues for at weak and strong 
QPC-qubit coupling respectively. The parameter <j> that pa- 
rameterises the scattering matrix equals n/16 (bottom) and 
4-7T-/5 (top) representing the weak and strong coupling lim- 
its respectively. For <f> = 7r/16 individual eigenvalues travel 
from 1 to cos7r/16 ~ 0.9808 at a rate of approximately 
one per 2it/eV. For <f> — 4tt/5, eigenvalues travel towards 
cos4-7r/5 ~ —0.8090 < at a rate of one per 2ir/eV, as shown 
in (a). Deviations from the correct asymptotics are due to 
finite size effects. Figure (b) contains the second derivative 

of P(t) = Det Q (0) (r) _1 4 (V) W- ( The second derivative is 
taken to remove an average slope and curvature.) Oscillations 
with period h/eV are seen (bottom) for cj> — tt/16, while for 
cj> — 4-7T-/5 (top), the periodicity of P(t) doubles. 



4 



minant P(r) is not precisely periodic in r since it decays 
at large r owing to decoherence. Still, it oscillates and 
the period of these oscillations doubles as we go from 
weak to strong coupling (Fig. 3b). The doubling can 
be understood in terms of the transfer of the eigenval- 
ues of Q( v '(t) upon increasing r (Fig 3a) assuming the 
parametrization ([5]). In the large r limit, energy-time un- 
certainty can be neglected in a "quasi-classical" approxi- 
mation: The operator II(t) projects onto a very long time 
interval, and is replaced by the identity operator. 
becomes diagonal in energy. All eigenvalues that are not 
equal to 1 are concentrated in the transport energy win- 
dow < E < eV where the filling factors in the QPC 
reservoirs are not the same. For s^si parametrized as 
in ((8]) these eigenvalues equal cos(^). There are eVr/2-K 
of them. In other words, the number of eigenvalues equal 
to cos cf> grows linearly with r. Numerical diagonalization 
of Q {v) (t) (Fig. 3a) shows that one eigenvalue is trans- 
fered from 1 to cos(4>) during time 2ir/eV. If cos (j> > as 
in the weak coupling case, this gives rise to P(t) oscilla- 
tions with frequency eV/2n manifesting integer charges. 
However cos <fi becomes negative at stronger couplings, 
so that P(t) changes sign with each eigenvalue transfer. 
Two eigenvalues have to transfer to give the same sign. 
The result is a period doubling of the oscillations in P(t) 
and hence half-integer charges. This resembles the be- 
havior of the wave vectors of the Bethe-Ansatz solution 

The parametrization ([8]) of the s 2 si is not general. 
However, the eigenvalue transfer arguments help to un- 
derstand general scattering matrices. Eigenvalue transfer 
still occurs at frequency eV/2Tr but instead of traveling 
along the real line, eigenvalues follow a trajectory inside 
the unit circle in the complex plane. Fractional charge is 
pronounced if the end point of the trajectory has a neg- 
ative real part. Numerical results for general scattering 
matrices are presented in Appendix B. 

Results presented so far are for "spinless" electrons. 
Spin degeneracy is removed by e.g. high magnetic field. 
If spin is included, but scattering remains spin indepen- 
dend, then two degenerate eigenvalues are transported 
simultaneously. In this case, the half-integer charge dis- 
sapears for the parametrization ([8|) but persists for the 
more general choice of complex eigenvalues. The results 
of further numerical work that confirm this are presented 
in Appendix C. 

We have studied a quantum transport setup that can 
easily be realized with current technology, namely that 
of a quantum point contact coupled to a charge qubit. 
The qubit is operated as a measuring device, its out- 
put signal — the probability p 2 — is directly seen in 
the QPC current. The dependence of the signal on the 
qubit level splitting reveals the nature of charged exci- 
tations in the voltage-driven QPC. When the qubit is 
weakly coupled to the QPC, the dependence reveals ex- 
citations with electron charge e. We demonstrated that 
for stronger coupling, the dependence suggests the exis- 
tence of the excitations that carry half the charge of an 



electron. 



APPENDIX A: NUMERICAL METHOD 

In this Appendix we give a more detailed account of 
the numerical calculation of the qubit tunneling rates 
ri2(e) and r2i(e) than is presented in the main text. 
Our starting point is Eq. (7) of the main text. In or- 
der to discuss qubit transitions from |1) to |2) as well 
as the reverse transition simultaneously, we change no- 
tation slightly. In what follows, indices i and / refer to 
the initial and final state of the qubit respectively. We 
consider "forward" transitions (/, i) = (2,1) and "back- 
ward" transitions (/, i) — (1,2). The central object of 
numerical work is the operator 



l + (sl§ f -l)iL(- T )f(V)( s ) T <0 



l + (4«i-l)n(r)/W(e) 



T > 



(Al) 



We recall that the matrices §i and s/ are the scattering 
matrices of QPC electrons when the qubit is in state i or 
/. n(r) is a time-interval operator, 



n(T) tMitV = S(t - i')<W | g 



0<t<T (A2) 
otherwise ^ ' 



f( v > (e) is diagonal in energy. It contains the filling fac- 
tors of QPC-electrons in the various channels, includ- 
ing any bias voltage that may be present. Its form 
in the time-basis (at zero temperature) is given be- 
low in Eq. (|A9|) . The operator QX\t) has an in- 
finite number of eigenvalues outside the neighborhood 
of 1 in the complex plain. This implies that a regu- 
larization of the determinant is needed. Indeed, if one 
naively assumes the unregularized determinant to be 
well-defined and possesing the usual properties of deter- 
minants, such as Det(AB) = Dct(A)Dct(P), one may 



show that 



Det Q<Y>(t) 



Were this true, 

it would have implied that T 12(e) = T2i(e). This cannot 
be correct. At low temperatures, the qubit is far more 
likely to emit energy than to absorb it, meaning that one 
of the two rates should dominate the other. 

Regularization is achieved by multiplying with the in- 
verse of the equilibrium operator. The operator Q/i(r) = 

Q%' 1 ( T )~ 1 Q%\ T ) on ly nas a finite number of eigenval- 
ues for finite r that are not in the neighborhood of 1, 
and so its determinant can be calculated numerically in 
a straight-forward manner. (In this expression, Q^) (r) is 

the operator Q when the QPC is initially in equilibrium, 
i.e. the bias voltage V is zero.) We therefore proceed as 
follows: We define 



P(r)=Det QgVr'CV) 



(A3) 



and P(e) = J dr e leT P{f) as its Fourier transform. The 
equilibrium rate T e ^(e) is known from the study of the 
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Fermi Edge singularity. It is 



we find a discretized operator 



r e f l(e) = h\ 2 e(-e fi ) w 



(A4) 



where E co is a cut-off energy of the order of Ep and 



1 

4^ 



Tr In 2 (s j Si) 



(A5) 



The logarithm is defined on the branch (— 7r,7r]. With 
the help of these definitions we have 



|^( £ ')P( £ - £ '), (A6) 



where our task is to calculate P(e) numerically 

The operator Q 2 W T ) will be considered in the time 
(i.e. Fourier transform of energy) basis. We restrict our- 
selves to the study of single channel QPC's, in which case 
the scattering matrices si and §2 are 2x2 matrices in 
QPC-channel space. We work in the standard channel 
space basis where 



rk A 
t k r' k 



with t, t' the left and right transmission amplitudes and 
r, r' the left and right reflection amplitudes. Because 
IT(t) is a projection operator that commutes with the 
scattering matrices, we can evaluate the determinant in 
the space of spinor functions tf)(t) defined on the interval 
t e [0, t]. (We consider r > 0.) Then 



Q^(t)V (*) 



where 



2S1- 



1) / dt'f^(t 



(A8) 



(s) 
8(eV-e) 

l-dv\ e- lteV - 1 



2n(t + i0^ 



2nt 



<A9) 



is the Fourier transform of the zero-temperature filling 
factors of the reservoirs connected to the QPC and + 
is an infinitesimal positive constant. Discretization of 
this operator proceeds as follows. We choose a timestep 
At < r such that N — t f At is a large integer. We 

will represent Q^if) (and Q^t) -1 ) as 2N x 27V di- 
mensional matrices. We define a dimensionless quantity 
r\ = eV At. P(t) can only depend on r in the combina- 
tion reV because there are no other time- or energy scales 
in the problem. We will therefore vary r by keeping N 
fixed and varying 77. Using the identity 



1 + (4si - i)n/ 

= S k i + (4*1 _ !) 



hi 



1 



+ ^( v 2^ 4 ' + 2tt»(Z - fc) {1 ~ Skl) 



(1 - <y«) 

(AH) 



nonequilibrium correction 

To test the quality of the discretization as well as its 
range of validity we do the following. When sjsi is close 
to identity, we can calculate P(t) perturbatively, both for 
the original continuous operators and for its discretized 
approximation. If we take s^si — e 
we find 



then to order 



Pcont.(r) = l + 2[j- 



N 



dz C0S{ ^ - 1 (N z) 



( Al2 ) 

where r = Nrj/eV for the continuous kernel while for the 
discretized version we find 



( A? ) ^disc» = l+ 2 



£) a g0OB(^-l (JV _. ) !A1>!) 



C=i 



C 2 



which indicates that the range of validity is 77 <g; 27r. 

In practice we take N = 2 s . Larger N would de- 
mand the diagonalization of matrices that are too large 
to handle numerically. We find results suitably accu- 
rate up to r\ = 7r/4, thereby giving us access to P(t) for 
|r| 6 [0,647r/eV]. 

To summerize, the procedure for calculating the tran- 
sition rates T21 and Ti2 is 

1. For given scattering matrices s\ and §2, calculate 
P{t) numerically using the discrete approximations 

for the operators Q 2 \\ T ) an d Q^i 1 ')- U se a fixed 
large matrix size, and work in units [r] = [ey] _1 . 
Generate data for many positive values of r. 

2. Extend the results to negative r by exploiting the 
symmetry P(t) = P(—t)*, and Fourier transform 
the data. 



3. Form the convolutions of Eq. IA6I with the known 
equilibrium rates to obtain the non-equilibrium 
rates. 



APPENDIX B: CHOICE OF SCATTERING 
MATRICES 

In the main text we confined our attention to the one 
parameter family of scattering matrices 



4*i 



cos 4 
i sin ( 



i sin < 

cos a 



(Bl) 
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For this choice, P(r) is a real function of time. For 
9 < tt/2 its fluctuations are associated with energies 
~ ±eV due to the transfer of eigenvalues from 1 to cos <j> 
at a rate of one per h/eV. For cf> > tt/2 however, cos <j> 
is negative and two eigenvalues have to be transfered be- 
fore the sign of P(r) returns to its initial value. The 
period of fluctuantions of P(r) doubles and becomes as- 
sociated with energies ±eV/2. Because P(r) is real, the 
fluctuations with positive and negative energies are equal: 
P(e) = P(— e). This translates into the following feature 
of the probability P2 to find the qubit in state |2). For 
4> < k/2, j>2 (e) changes from 1 to in an energy inter- 
val of length 2eV. For tj> > tt/2, this interval shrinks to 
eV. The boundry of the interval is defined more sharply 
the closer is to or tt. The shrinking from 2eV to 
eV of the interval in which p2 varies significantly is ex- 
plained in terms of charge fractionalization: For <j) > tt/2 
the excitations in the QPC transmit half the charge of 
an electron so that the energy that the qubit can absorb 
from the QPC changes from eV to eV/2. 

Since the QPC scattering matrices contain parameters 
that are not under experimental control, it is relevant 
to ask how the results are altered when a more general 
choice 



4*1 = 



e cos < 
i sin <b 



e %e cos < 



(B2) 



with 4> <E [—§,§] and 9 £ [0,tt] is made for the scatter- 
ing matrices. With this choice, eigenvalues travel from 
I to e cos</> at a rate of one per 2n/eV. This means 
that the period doubling of P(r) no longer takes place. 



The phase of P(r) does not return to its original value 
after the transfer of two eigenvalues. Rather, one expects 
fluctuations associated with an energy (n — ^)eV, n = 
0, ±1, ±2, . . . Because P(r) is no longer real, positive and 
negative frequencies don't contribute equally. However, 
while the eigenvalue trajectories lie close to the real line, 
one can expect results similar to those obtained for real 
P(t). We obtained numerical results for four scattering 
matrices of the form (IB2 



We chose 9 = ^tt, 



1 7r, §7r and 
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|tt. To sharpen abrupt features we chose <f> = tt/9 so 
that the exponential decay of P(r) is associated with a 
long decoherence time: ~ 0.06U/eV. As depicted in Fig. 
0, we found P(e) to behave as follows. For 9 close to 



zero, P(e) consists of one peak situated at e 



2tt 



eV. 



The tails of this peak vanish at e — (±1 — eV. The 
closer to zero that 9 is taken, the more abrupt this be- 
havior of the tails become. As 9 is increased, a sec- 
ond peak starts appearing at e = (f — 4A eV. When 
9 = tt — + , the height (and width) of this peak ex- 
actly equals that of the peak at — -^eV. In the interval 
e E [—-^eV, (f — j-) eV] that is bounded by the peaks, 
P(r) is significantly larger than in the region outside the 
peaks. This behavior of P(e) translates into the occupa- 
tion probabilities p2(e) depicted in Fig. ((5]). For 9 < tt/2, 
P2{£) still changes from unity to zero in an interval of 
length 2eV manifesting excitations with charge e while 
for 8 > tt/2 the interval shrinks to eV, indication half- 
integer charge. The closer 9 moves to or tt, the sharper 
the interval becomes defined. We therefore conclude that 
the fractional charge phenomenon in the QPC is not con- 
fined to the special choice (|B1|) of scattering matrices. 
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Level splitting, e [eV] 
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FIG. 4: The function P(e) that contains the effect of 
the bias voltage V. As explained in the text, sjsi was pa- 
rameterized as in Eq. (|B2|l . A value <f> — ? is used througout. 
The values of 9 in (a), (b), (c) and (d) are respectively f , f, 
^ and When 6 < tt/2, then P(e) has a fairly symmetric 
peak centered at —eV6/2n. The tails of this peak vanish at 
e ~ (—8/2n ± l)eV. When 9 > tt/2, there are two asymmet- 
ric peaks at -eV8/2iT and (I - 8/2ir)eV. The value of P(e) 
is significantly larger for e G [— eV8/2ir, (I — 8/2Tr)eV] than 
outside this interval. 



FIG. 5: The probability p-i{s)- s\s\ is chosen as in Fig. 
A value 4> — f is used througout. The values of 9 in (a), (b), 
(c) and (d) are respectively f , f , ^ and ^f. When 8 < tt/2, 
the occupation probability p2 is significantly different from 
its asymptotic values and I in an e interval of 2eV. When 
8 > tt/2, this interval shrinks to eV. The boundaries of the 
interval are more sharply defined the closer 8 is to n/2. The 
shrinking of the interval corresponds to a cross-over in the 
QPC from excitations that transmit charge e to excitations 
that transmit charge e/2. 
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APPENDIX C: INCLUSION OF SPIN 

Up to this point we have considered spinless electrons 
in the QPC. In this Appendix we investigate the effect of 
including spin. We still take the interaction between the 
QPC and the qubit to be spin independent. However, 
the mere existence of a spin degree of freedom for QPC 
electrons doubles the dimension of channel space. The 
narrowest QPC now has two channels in stead of one and 
P_i(r) = -P s= o(t) 2 , i.e. the determinant P _i(t) with 

spin included is the square of the determinant P s= o(r) 
without spin. For real determinants, squaring kills the 
phase. This means that the observed period doubling 
for the parametrization of Eq. (|B1[) disappears and with 
it the half integer charge features of P2- Physically, it 
could be that two charge e/2 excitations are transmitted 
through the QPC simultaneously. However, fractional 
charge is saved by the fact that, for 9 / 0, P s= o(e) has 
two peaks with different heights. Suppose the relative 
peak heights are A and 1 — A, i.e. 

P s=0 (r) - (1 - A)e l ^ eVT + Ae- {1 -^ )eVT (CI) 

where A is a real number between and \. (A = 
corresponds to 9 = while A=\ corresponds to 9 = tt.) 
It follows that P _i (s) has three peaks at 

1. e = -2^eV with height (1 - A) 2 , 

2. e = (l - 2^) eV with height 2A(1 - A) 



3. and e = (2 - 2^) eV with height A 2 

As long as A is small, i.e. 9 is not too close to tt, the 
first two peaks will dominate the third, and a signature 
of fractional charge may still be observable in P2(s)- Fig. 




- 1 1 - 1 1 




- 1 1 - 1 1 

Level splitting, s [eV] 



FIG. 6: The probability P2(e) with spin included, s^si is 
chosen as in Fig. Q and ([5]): A value (f> — ^ is used througout. 
The values of 9 in (a), (b), (c) and (d) are respectively f , f, 
^j- and Fractional charge features are still clearly visible 
for 9 > tt/2. 

(|6|), contains P2 calculated for the same scattering ma- 
trices as in Fig. (|5|), but with spin included. The cases 
when 9 — %tt and 9 = |7r still contain clear half-integer 
charge features. For 9 very close to tt (not shown) these 
features disappear. 
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